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Abstract 

We investigate the pp — > pprj and pn — * pnr\ reactions within an effective Lagrangian model 
for laboratory kinetic energies ranging from very close to the eta meson production threshold to 
about 10 GeV. Production amplitudes include contributions from the mechanism of excitation, 
propagation and decay of iV*(1535), A r *(1650), and A?"* (1710) baryonic resonances. The initial 
interaction between two incoming nucleons is modeled by the exchange of tt, p, u and a mesons 
where the vertex parameters are taken to be the same as those used in the previous applications 
of this model. Parameters of the resonance vertices also have been taken from our earlier studies 
wherever applicable. Calculations have been done for total as well as differential w production 
cross sections. To describe the data for energies closer to the production threshold final state 
interactions among the out-going particles have been included by means of a generalized Watson- 
Migdal method. Terms corresponding to the excitation of iV*(1535) resonance and pion exchange 
process dominate the cross sections. With a single set of vertex parameters our model is able to 
describe the available data well on total cross sections for beam energies ranging from close to 
threshold to upto 10 GeV. 

PACS numbers: 13.60.Le, 13. 75. Cs, 11.80. - to, 12.40. Vv 



I. INTRODUCTION 



The low energy behavior of the quantum chromodynamics (QCD) is not accessible to 



the p— e approaches; the Lattice g ao g e theo ry Q is the idea, too, tor this puIp ose. 
Despite enormous computational power necessary for the numerical realization, lattice QCD 
calculations have started, very recently, to describe masses and other constants of the bary- 
onic ground as well as excited states [2J. Experimentally, the determination of baryonic 
resonance properties proceeds indirectly by exciting the nucleon with the help of a hadronic 
or electromagnetic probe and performing measurements for their decay products (mesons 
and nucleons). The reliable extraction of nucleon resonance properties from such experi- 
ments is a major issue of the hadron physics. 

In recent years, important advances have been made in the experimental investigation 



of meson production reactions in nucleon-nucleon (NN) collisions, particu 



energies very close to respective production thresholds 
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151 ] . Low incident energies also provide the opportunity to investigate the meson- nucleon 



interactions through these reactions since in this energy regime the final state interactions 
among the outgoing particles affect strongly the meson production cross sections. 

The r) meson, which is the next lightest nonstrange member in the meson mass spectrum, 
has been a subject of considerable interest. It has been thought of as a probe for the ss 
component in the nucleon wave function [0]. There is also interest in measuring the rare 
decays of rj which could provide a new rigorous test of the standard model 13] or even of 
the physics beyond this. The nucleon resonance iV*(1535) [S , n(1535)] with spin ~, isospin 
|, and odd parity, has a remarkably large r/N branching ratio. It lies very close to the 
threshold of the NN —>■ NNrj reaction and contributes to the amplitude of this reaction 
even at the threshold. Therefore, the study of //-meson production in NN collisions at the 
near threshold beam energies provides the unique opportunity to investigate the properties 
of iV*(1535) which have been the subject of some debate recently (see, e.g., [l8|] ) . The 
attractive nature of the 7/-nucleon interaction may lead to the formation of bound (quasi- 
bound) 77-nucleus states (see, e.g., l2lL I22L |23j]). This subject has been a topic of 



intense discussion at a recent work shop 



24]. 



Production of rj meson in heavy ion collisions also is of great interest. Due to the high 
threshold of the elementary production reaction, rj mesons in such collisions are produced 
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only by very energetic nucleons and reflect therefore, the tails of the nucleon momentum 



distributions as they arise in a high density and high temperature phase of the collision [25]. 
The elementary NNrj production cross sections are a crucial ingredient in the transport 
model studies of the ^-meson production in the nucleus-nucleus collisions. 

Since lattice QCD calculations are still far from being amenable to the low and interme- 
diate energy scattering and reaction processes, the effective methods are mostly used for the 
description of the dynamics of the meson production reactions in hadronic collisions. These 
approaches introduce the baryonic resonance states explicitly in their framework and try to 
extract their properties by comparing the theoretical results with experimental observables. 



Several authors have used mode 
collisions 
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27, 



28 



29 



30 



31 



s of such type in describing the r] meson production in NN 



32]. 



The main objective of this paper is to investigate the r] meson production in NN collisions 
in the framework of an effective Lagrangian Model (ELM) which has been used earlier rather 



successfully to describe the pion [33|, |34j, associate kaon [351], and dilepton [36| production 



data in such collisions. The motivation here is to see as to how far can one exp 



recently measured data on total and differential cross sections of pp — > ppt] 
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ain the 

,y,Q 

and the pn — > pnrj [l0[ reactions within this model using same sets of entrance channel and 
resonance channel parameters (for those resonances that appeared in earlier applications 



35, 



36|). 



Within the ELM, initial interaction between two incoming nucleons is modeled by an 
effective Lagrangian which is based on the exchange of ir, p, to and a mesons. The coupling 
constants at the nucleon-nucleon-meson vertices are determined by directly fitting the T 
matrices of the NN scattering in the relevant energy region 37j. The effective Lagrangian 



uses the pseudovector (PV) coupling for the nucleon-nucleon-pion vertex, as it is consistent 



with the chiral symmetry requirement of the QCD [38|, [39|, |40| and also it leads to negligible 
contributions from the negative energy states ("pair suppression phenomena") [4l[. The rj 
meson production proceeds via excitation of iV*(1535), N*(1650) and iV*(1710) intermediate 
baryonic resonance states which have known branching ratios for the decay into the r]N 
channel. The coupling constants at the resonance-nucleon meson vertices are determined 
from the experimental widths for the decay of the resonances into the relevant channels 
except for those involving the to meson where they are determined from the vector meson 
dominance (VMD) hypothesis. The interference terms between amplitudes corresponding 
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to various meson exchanges and the intermediate resonance states have been included. 



The final state interaction (FSI) among the outgoing particles affect str ong ^ 
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y the cross 



321 ] . In ap- 



sections of the NN — > NNr] reaction at near threshold beam energies 
plications jjf the ELM to describe the near threshold meson production reactions in NN 
collisions 



34 



351 ]. the FSI effects were included within the Watson-Midgal theory [42] which 



is based on the assumption that the FSI effects are strong in relation to the production 
process and that they occur attractively between only one particular pair of out-going par- 
ticles. In Ref. 35j, this method was used somehow arbitrarily for all the 3 outgoing pairs 



of particles. However, Watson's method as such is not applicable, in strict sense, for those 
cases where the attraction between outgoing particles is not so pronounced or where in- 
teraction between more than one pair is to be included in calculations. In this paper, we 
employ a generalized Watson method in which three-body states are treated by splitting 
the total interaction into pairwise net interactions which leads to a series decomposition of 
the net scattering among all the particles in terms of separate total scattering between pairs 
of particles (see, e.g. 43]). However, the three-body interactions are neglected. In view of 



the arguments presented in Refs. 32|, |44| in favor of using the three-body scattering theory 



to describe the NNr] process, it would be interesting to see to what extent this generalized 
method is able to explain the NNr] production data. 

In the next section, we present a brief description of the ELM where we describe the main 
ingredients of the theory and give all the input parameters used in our calculations. The 
generalized Watson method of FSI effects is also described in a subsection here. The results 
of calculations and their discussions are presented in section III. Summary and conclusions 
of our work are given in section IV. 



II. FORMALISM 



A representative of the lowest order Feynman diagrams contributing to the r] meson 
production considered by us, is shown in Fig. 1. Momenta of various particles are indicated 
in Fig. la. q, Pi, and p v are four momenta of the exchanged meson, the intermediate 
resonance and the r^-meson, respectively. To evaluate various amplitudes, we have used the 
effective Lagrangians for the nucleon-nucleon-meson and resonance-nucleon-meson vertices 
as described below. 
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FIG. 1: Feynman diagrams for the ?7-meson production in nucleon-nucleon collisions. Diagrams 
(a) and (b) show the direct processes, while (c) and (d) the exchange ones. R represents a baryonic 
resonance. 



Nucleon-nucleon-meson vertex 



As done before in the investigation of pp — > ppir , pp — > pmr + [33], pp — ► pAK + 35], 
and NN — > NNe + e~ 36j| reactions, the parameters for nucleon-nucleon-meson vertices are 
determined by fitting the iViV elastic scattering T matrix with an effective NN interaction 
based on the jr, p, u, and a meson exchanges. The effective NNmeson Lagrangians are 



v see, e.g. 
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C>NN-k 
C-NNp 

£>NNw 
£-NNa 



9Nm 
2m 



N 



-9NNp^N hp + ^^ a ^ 9 j T ■ P^X- 



-gNNw^N yip + 

Qnn^nct^n- 



2m N ) 



(1) 
(2) 

(3) 
(4) 



In Eqs. (l)-(4), we have used the notations and conventions of Bjorken and Drell j47| and 
definitions of various terms are the same as those given there. In Eq. (1) denotes the 



nucleoli mass. It should be noted that we use a PV coupling for the NNtt vertex. Since 
these Lagrangians are used to directly model the NN T matrix, we have also included a 
nucleon-nucleon-axial-vector-isovector vertex, with the effective Lagrangian given by 

C-NNA = 9NNA^l5l^ ■ A**, (5) 

where A represents the axial-vector meson field. This term is introduced because if the mass 
of the axial meson A is chosen to be very large tun) [stJ and gNNA is defined as 



1 / gNNn\ ,„x 
—f^A 7: > ( 6 ) 



9nna = ^) 

it cures the unphysical behavior in the angular distribution of NN scattering caused by 
the contact term in the one-pion exchange amplitude. It should be mentioned here that A 
is different from the ai(1260) meson resonance. The role of the A vertex is to explicitly 
subtract out the contact term of the one-pion exchange part of the NN interaction. Similar 
term in the coordinate space potential is effectively switched off by the repulsive hard core. 
At each interaction vertex, the following form factor is introduced 



,-iAW 



,i = ir,p,a,u}, (7) 



1 \ \ 2 a 2 / ' 

where and are the four momentum and mass of the ith exchanged meson and Aj is 
the corresponding cut-off parameter. The latter governs the range of suppression of the 
contributions of high momenta which is done via the form factor. Since NN scattering cross 
sections decrease gradually with the beam energy (beyond certain value), and since we fit 
the elastic T matrix directly, the coupling constants are expected to be energy dependent. 
Therefore, while keeping the cut-offs A, [in Eq. (7)] energy independent, we take energy 
dependent meson-nucleon coupling constants of the following form 

9(\^) = 9oexp(-£y/s), (8) 

where s is the square of the total CM energy. The parameters go, £, and A were determined by 
fitting to the T matrix of the relevant proton-proton and proton-neutron scattering data at 
the beam energies in the range of 800 MeV to 4.00 GeV 37j]. This procedure also fixes signs 
of the effective Lagrangians [Eqs. (l)-(5)]. The values of various parameters are shown in 
Table I [the signs of all the coupling constants (g) are positive]. In this table the parameters 
of the A exchange vertex are not explicitly shown as they are related to those of the pion 



via Eq. (6). We would like to remark that the same parameters were also used to describe 
the initial NN interaction in the calculations reported in Refs. 33j, |35|, |36j. This ensures 
that the elastic NN elastic scattering channel remains the same in the description of various 



TABLE I: Coupling constants for the AA-meson vertices used in the calculations 



Meson 


g 2 /An 


I 


A 


mass 








(GeV) 


(GeV) 


7T 


12.562 


0.1133 


1.005 


0.138 


a 


2.340 


0.1070 


1.952 


0.550 




46.035 


0.0985 


0.984 


0.783 


P 


0.317 


0.1800 


1.607 


0.770 


k p = 6.033, k u = 0.0 











inelastic processes within this model. 

The main criterion for choosing the meson exchanges as discussed above is to describe the 
NN scattering in the relevant beam energy region. We have left out the ?7-meson exchange 
in our description of the NN interaction since due to the psuedoscalar (PS) nature of its 
coupling, the contribution of the exchange terms of this particle (having a mass much larger 
than that of the psuedoscalar pion) is ex pec ted to be very small. Furthermore, as confirmed 



by several studies (see, e.g., 
rather small. 



48 



49 



50 



511]). the coupling constant for the NNr] vertex is 



B. Resonance-nucleon-meson vertex 



As the rj meson has zero isospin, only isospin-| nucleon resonances are allowed. Below 
2 GeV center of mass (cm.) energy, AT* (1535) has a prominent decay branchin g ra tio of 
40-60% into the Nr] channel (see, e.g., the latest review of the particle data group |52j). On 
the other hand, A/"* (1650) and AT*(1710) resonances have a small but non-negligible decay 
branching ratios of 3 - 10% and 6 ± 1%, respectively, to this channel. In comparison to 
these, the branching ratio for the decay of A^*(1520) resonance to Nrj channel is negligibly 
small and we have not included it in out description. In several previous studies of the 
AW — > NNr] reaction, contributions from only AT* (1535) resonance have been considered. 



Since all the three resonances can couple to the meson-nucleon channels considered in 
the previous section, we require the effective Lagrangians for all the four resonance-nucleon- 
meson vertices corresponding to all the included resonances. At the spin-| resonances - 
N — 7r (77) vertices, we have the choice of RS or PV couplings. The corresponding effective 
Lagrangians can be written as 35, 



48 



53, 



54| 



r PV 
r PS 

rPV 
r PS 



9N* Nit - 

1/2 (0^)^ + 11.(5. 



M 

■gNi /2 Nn^N*iFr® n ^ N + h.c, 

9N* n _ 

' "^f^T-^^^ + h.C, 



M 

-gN*„Nr,^N*iFT& v ty N + h.C, 

1/2 



(9) 
(10) 

(11) 
(12) 



where M = (m^* ± mjv), with upper sign for even parity and lower sign for odd parity 
resonance. The operators r, r M , are given by, 

r = 75, r M = 757^, (13) 
r = 1, r M = 7m , (14) 

for resonances of even and odd parities, respectively. We have performed calculations with 
both of these couplings. The effective Lagrangians for the coupling of resonances to other 
mesons are, 

1 



Cm* 



5 *jv.- T^d u r ■ p^ N . + h.c. 



■ /2 N P ~ -g N1/2 N p VN*—L^ 

-9n* ,n^n*~ r^d u uj^ N . + h.c. 

£n* /2 n<t = 9N* /2 N^N*^'cr^N + h.c, 



C N ; /a Nu 



where operators V and r„„ are 



r' 



1 r 



a 



fJ,V 



r' = 75, r 



[XV 



75<v> 



(15) 

(16) 
(17) 

(18) 
(19) 



for resonances of even and odd parities, respectively. 

We assume that the off-shell dependence of the NN* vertices are determined solely by 



multiplying the vertex constants by form factors. Similar to Refs. 
following form factors for iV*iV-nieson vertices 

(Af) 4 



46 



55(, we use the 



(Af *) 4 + (qf - mf ) 2 



(20) 
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The resonance couplings are determined from the experimentally observed branching 
ratios for the decay of the resonances to the corresponding channels. Since the resonances 
considered in this study have no known branching ratios for the decay into the Nuj channel, 
we determine the coupling constants for the N*Nu vertices by the strict vector meson 
dominance (VMD) hypothesis [56|], which is based essentially on the assumption that the 
coupling of photons to hadrons takes place through a vector meson. For details of these 
calculations we refer to j35|. 

The resonance properties and the values of various coupling constants are given in Table 
II. Value of the cut-off parameter (Aj^* ) is taken to be 1.2 GeV for all the vertices, which is 



36l . |55| |. Fixing A^ to one value minimizes the number of 



the same as that used in Refs. 
free parameters. 

It should however, be stresses that the branching ratios determine only the square of 
the corresponding coupling constants, thus their signs remain uncertain in this method. 
Predictions from independent studies are used as a guide to fix these signs. We have followed 
here the results of Ref. 4(| for this purpose. The propagators for various mesons and nucleon 
resonances in the calculation of the amplitudes have been taken to be the same as those 



discussed in 



35 



36|. 



C. Amplitudes and cross sections 

After having established the effective Lagrangians, coupling constants and form of the 
propagators, the amplitudes for various diagrams associated with the NN — » NNrj reaction 
can be calculated in straight forward manner by following the well known Feynman rules. 
The isospin part is treated separately. This gives rise to a constant factor for each graph, 
which is shown in Table III. 

It should be stressed here that signs of various amplitudes are fixed, by those of the 
effective Lagrangians, coupling constants and propagators as described above. These signs 
are not allowed to change anywhere in the calculations. 

The general formula for the invariant cross section of the NN — > NNrj reaction is written 

as 

^ = nr, m% , T^ ^( p f ~ Pi)W A % (2D 



2j[( Pl - V 2?-m%} W L\ E a 
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TABLE II: Resonance parameters and the coupling constants for various decay vertices. Coupling 
constants at the N*Nco vertices are obtained from the vector meson dominance hypothesis (see, 
e.g. Ref. [35j). 



Resonance 


Width 
(GeV) 


Decay channel 


9 


iV*(1535) 


0.150 


iVvr 


0.6840 






Np 


3.9497 






Nto 


1.4542 






Na 


2.5032 






Nrj 


2.2000 


iV*(1650) 


0.150 


Nir 


0.8096 






Np 


2.6163 






Nlo 


1.8013 






Na 


2.5032 






Nr) 


-0.5469 


iV*(1710) 


0.150 


Nir 


1.0414 






Np 


2.9343 






Nlo 


1.5613 






Na 


0.6737 






Nr] 


1.0328 



where Tji represents the total amplitude, Pi and Pf the sum of all the momenta in the initial 
and final states, respectively, and p a the momenta of the three particles in the final state. 
The corresponding cross sections in the laboratory or center of mass systems can be written 
from this equation by imposing the relevant conditions. 

D. Final state interaction 

For describing the data for the NN — > NNr] reaction at beam energies very close to 
the j] production threshold, consideration of the final state interaction (FSI) among the 
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TABLE III: Isospin factors for various diagrams, Isovector corresponds to it and p exchange graphs 
while isoscalar to uj and a ones 



graph 


isovector 


isoscalar 




TtTI — > T)T)Ti 




direct 


1.0 


1.0 


exchange 


1.0 


1.0 




pn — ► pnr\ 




direct 


-1 


1 


exchange 


2 






three out going particles is important. We follow here an approximate scheme in line with 
the Watson-Migdal theory of FSI 42| . In this approach the energy dependence of the cross 
section due to FSI is separated from that of the primary production amplitude. This is based 
on the assumption that the reaction takes place over a small region of space, a condition 
fulfilled rather well in near threshold reactions involving heavy mesons. This method has 
been extensively applied to study the low momentum behavior of the pion 34], |57|, |58j, r\ 



meson 



59 



60 



611 ]. associated hyperon 35|, |62) and 0-meson 63j production in NN collisions. 



The total amplitude is written as 

T fi = T (NN — > NNrj) ■ Tff, 



(22) 



where T (NN —>■ NNrj) is the primary production amplitude, while Tff describes the re- 
scattering among the final particles which goes to unity in the limit of no FSI. The factoriza- 
tion of the total amplitude into those of the FSI and primary production (Eq. (22)), enables 
one to pursue the diagrammatic approach for the latter within an effective Lagrangian model 
and investigate the role of various meson exchanges and resonances in describing the reac- 
tion. 



Watson's original method 42j was developed for those cases where the final state inter- 
action is strong in relation to the production process and where it is confined only to one 
particular pair of particles (mostly among nucleons in case of nucleon-nucleon-meson final 
states). On the other hand, in certain cases it may be necessary to include FSI among all 
the three outgoing particles since even if the meson-baryon interactions are weak, they can 
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still be influential through interference. In Ref. [35j, the T matrix Tff was written (without 
providing any proof) as a coherent sum of the transition matrices describing the final state 
interaction among all three two-body subsystems of the final state nucleon-nucleon-meson 
system. We show here that this result can be obtained (in a slightly different form) by 
following the technique of multiple FSI as discussed in Ref. [43]. To that end, we first try 
to get a representation for the total amplitude T/j in terms of an expression similar to that 
given by Eq. (22) where the FSI amplitude Tff is appropriately constructed. 

To introduce the treatment of the FSI for the three-particle system, we assume that 
the three-particle final states can be represented by additive potentials of the form U = 
U12 + U 3 i + U23 = U 3 + Ui + U\. With this assumption, the total amplitude T/j can be 
written as a iteration series in terms of the production amplitude T (defined in Eq. (22)) 
and the three-body final state pair interaction amplitudes Tj = (see Appendix A for 
details). 

Tfi = To + J^TkGoTo + J^TkGoTjGoTo + ...., (23) 
k k^j 

where the final state interaction transition matrices Tk are as defined in Appendix A. Go 
is the Green's function corresponding to the free Hamiltonian (kinetic energy). Neglecting 
processes depicted by the third term, this decomposition can be expressed by Fig. 2. It is 
easy to show that the result of the Watson FSI theory are recovered if one retains only the 
amplitude T and a single pair amplitude, say, Ti 2 . 

Any practical calculation requires evaluation of the matrix element < 
123\Ti2GqTq\NiN2 >, where Ni and N2 denote two particles of the incident channel 
and 1,2 and 3 represent the outgoing channel particles. Introducing a complete set of 
intermediate states of particles, say, 1 and 2, we get for this channel 

1 rd 3 k[d s k' 3 N <k 1 k 2 |T 12 |k' 1 k' 2 xk' 1 k^k 3 |T |Ar 1 Ar 2 > 
L tt = T^Tv? / -WWTWET 6 l k i + k 2 - k i - k 2 



(2tt) 3 J 2E[ 2E' 2 v 1 1 zy E — (E 3 + E[ + E' 2 ) + ie 

(24) 

It may be noted that in this form, each particle lies on the mass shell and three momentum 
is conserved in the intermediate processes. Introducing the total and relative momenta, 
p' = k' x + k' 2 , 2q' = k' x — k 2 , and evaluating the integral in the barycentric frame of 1 and 2, 
we obtain, 

„ 1 f r 12 (£, 6- 6')T (?, 9', k 3 ; kj 2g'#W 

T ^-WfJ T^TTe ^— • (25) 

12 



FIG. 2: Final state scattering among three particles evaluated to the lowest order 

where £ is the energy of 1 and 2 in this frame, £' is the intermediate state and denotes the 
orientation (9, <p) of q with respect to a fixed axis. 

For further evaluation of the integral, we make a partial wave decomposition of the 
amplitude and for each partial wave we rewrite it as 



Now we make the assumption that the ratio within the square brackets in the integrand of 
Eq. (26) is constant upto a certain energy £ c and zero thereafter. Extending this procedure 
to all the three interacting pairs, we get in the low energy and s-wave limit 



TfT oc T 12 (£,OT (0 / 



[Tufa Qro(Q/ri2(g, Qr (fl] 2 q >d? ( m> 



(26) 



(27) 



In Eq. (27) T ff is defined as 



Tff(0 = Y.'uDiiCa- 



(28) 



where 




(29) 
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It is obvious that Eq. (28) allows interference among the final state scattering amplitudes. 
We further note that apart from the factor c^, this equation is similar to that used in 



Ref. 



35]. 



The derivation of Eq. (27) is independent of the strength of the interaction and of whether 
it is attractive or repulsive. The quantity c^- can be regarded as the amount of final state 
scattering that takes place in a particular channel ij. A plausible value of the cut off £ c 
comes from the constraint that for the NN sub-state should come out to be unity since 
in the limit of FSI in only this sub-state, we should recover Watson's result. With the same 
value of the £ c , the Cjj for the r/N sub-state comes out to be 1.07. It must however, be noted 
that this procedure determines the value of the c^- at best only for the NN channel. It 
remains largely undetermined for the rj — p sub-state which could even be dependent on the 
relative energy of this channel. Since, for the time being we do not have a definite method 
to determine this quantity, we have taken the same value for the parameter £ c for both the 
sub-states. 

For calculating the FSI amplitude for the rjN sub-state, we note that there are no direct 
measurements of the elastic r]N scattering and the information about the rjN elastic scat- 
tering amplitude is obtained by describing the nN — > rjN and 7iV — * rjN data within some 
model. Recently, there has been suggestions to determine the rjN scattering amplitude from 
the studies of associated photoproduction of (ft and rj mesons off the proton [64j. We adopt 



here the results reported in Ref. 65[ where rjN scattering parameters have been obtained 
by fitting the irN — > ttN, ttN — > rjN, ^N — > rjN data in an energy range from threshold to 
about 100 MeV, in a K matrix method. These authors write the elastic rjN scattering T 
matrix as 

T~ l = l/a + r -^q 2 v + sq* - iq v , (30) 
where q v is the momentum in the rjN center of mass (cm.) system. Seven sets of values 



of the parameters a, r and s are given in Ref. [65|]. We found that the best description 
of the data (within the realm of our overall input parameter sets given in Tables I and 
II), is provided by the r]N scattering amplitudes with the parameter set, a = 0.51 + i0.26 
fm, ro = —2.50 — z0.310 fm, and s = —0.20 — i0.04 fm 3 . It should be noted that the real 
part of the rjN scattering length of this parameter set (or) is about half of that of the 
"preferred" set of Ref. [65|. A larger a# is also supported by the calculations presented in 
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FIG. 3: The total cross section for the p + p — > p + p + r/ reaction as a function of the beam 
momentum. The dashed, dotted and dashed-dotted curves represent the contributions of N* (1535), 
iV*(1650), iV*(1710) baryonic resonance intermediate states, respectively. Their coherent sum is 
shown by the solid line. The experimental data are from [741 ] . 



Ref. 



661 ] . However, we note that in the theoretical description of the pp — > pprj reaction as 



reported in Refs 
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671 ] . the value of q r w as similar to that used by us. A smaller clr is also 



consistent with that extracted in Ref. 



54J ] in an effective Lagrangian model analysis of the 



meson-nucleon scattering. Furthermore, it was noted in Ref. 68( that within a three-body 
model, the shapes of the np — > rjd cross sections can be explained over a wide energy range 
only with a cir around 0.42 fm. A smaller value of Or is also consistent with the Jiilich 



model 



69] 



The FSI amplitude Tnn has been calculated by following the Jost function method using 
the effective range expansion (ERE) of the iViV phase shifts as discussed in Refs. 34j.l35l.l42]. 
In case of the proton-proton sub-state, the Coulomb modified effective range expansion has 
been used 70|]. The effective range parameters for the iViV channel have been taken to be 



the same as those used in Ref. I3j 



It should however, be mentioned here that the use of on-shell forms to describe the FSI 
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T matrices Tit has been criticized by some authors. It has been argued in Refs. 7l|, [72 



that the absolute magnitudes of the cross sections obtained by such a procedure could be 

n 

uncertain because of the off shell effects. Even the Jost function method has been shown [73J 
to produce inadequate results in an study where the scattering length parameters for the 
A — p final state has been extracted from the pp — > pAK + data. In the next section we have 
examined the role of the off shell effects in the A A sub-state in some more details. 



III. RESULTS AND DISCUSSIONS 

The major aim of this paper is to check the suitability of our model and the vertex 
parameters appearing therein to describe the r\ production cross sections over a wide range 
of beam energies. We have therefore, applied our approach to describe the total cross sections 
for the pp — > ppr) reaction for beam energies ranging from near threshold to upto 10 GeV, 
and for the pn — > pnr) reaction for beam energies from threshold to upto 1.6 GeV. These 
are the energy regimes in which experimental data are available for the two reactions. We 
have also used this method to describe one set of the exclusive data, namely the r\ angular 
distributions for the former reaction. Calculations have been performed by using both the 
PS and PV couplings for N*Ntt and N*Nr] vertices. We note that the cross sections remain 
almost unchanged by switching from one type of coupling to another. In all the calculations 
shown below the coupling constants and cut off parameters for various vertices were the 
same as those discussed in section II. 

A cleaner check of the vertex parameters used in calculating the amplitude T (NN — > 
NNrj), is provided by the comparison of our calculations with the data for beam momenta 
above 3 GeV/c, since at these energies FSI effects are most likely to be unimportant. In 
Fig. 3, we show the comparison of our calculations and the experimental data (taken from 
Ref. (74|) for the total cross section of the pp — > ppr/ reaction at higher beam energies. We 
notice that the measured cross sections are reproduced reasonably well by our calculations 
(solid line) in the entire range of beam momenta. 

Individual contributions of various nucleon resonance intermediate states to the pp — > ppr] 
reaction are also shown in Fig. 3. Cross sections corresponding to A* (1535), A* (1650) and 
A* (1710) resonances are represented by dashed, dotted and dashed-dotted lines, respectively 
while their coherent sum is shown by the solid line. We note that the contributions of 
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FIG. 4: The total cross section for the p+p — ► p+p + ij reaction as a function of the excess energy. 
The dotted and dashed curves represent cross section obtained with FSI effects included only in the 
proton-proton sub-state of the final channel and no FSI at all, respectively. The solid line shows 
the results obtained with full FSI effects included as discussed in section II. The experimental data 
are from ll|, |12|, \l 



the iV*(1535) resonance dominate the total cross section for all the beam momenta. In 
comparison, those of iV*(1650) and AT* (1710) resonances are smaller by factors ranging from 
5-10. However, the interference terms of the amplitudes corresponding to various resonances 
are not negligible. It must again be emphasized that we have no freedom of choosing the 
relative signs of the interference terms. 

The results shown in Fig. 3 fix the parameters of all the vertices. In the application of 
our model to describe NNr] data at near threshold beam energies, the amplitude T (NN — > 
NNrj) has been calculated with exactly the same values for all the parameters. For these 
energies the FSI effects in the outgoing channels, have been included by using Eqs. (27-30). 
The experimental cross sections in this energy regime are given as a function of the excess 
energy (e) which is denned as e = y/s — 2m p — m v , where \fs is the invariant mass. 

In Figs. 4 and 5, we present comparisons of our calculations with the experimental data for 
total cross sections of the pp — > ppr\ and pn — ► pnrj reactions, respectively, as a function of e. 
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FIG. 5: The total cross section for the p + n — > p+n + n reaction as a function of the excess energy. 
The dotted and dashed curves represent cross section obtained with FSI effects included only in the 
proton-proton sub-state of the final channel and no FSI at all, respectively. The solid line shows 
the results obtained with full FSI effects included as discussed in section II. The experimental data 
are from 10[ ] . 



The solid, dashed and dotted lines show the results obtained by including the full FSI effects 
in all the three sub-systems, FSI only in pp and pn channel and no FSI at all, respectively. It 
should be noted that no arbitrary normalization constant has been introduced in any of the 
results shown in these figures. For the case of the pp — > ppn reaction, our full calculations 
describe the data quite well for e values in the range of 15 - 130 MeV. On the other hand, 
for the pn —> pnrj reaction they are in excellent agreement with the available data for all 
the beam energies. The FSI in the r]p sub-state is indeed quite important in our model. 
The difference between results obtained with FSI in only the pp sub-state and that in all 
the three subsystems of the final channel, is comparable to that reported in the three body 
calculations of Ref. These authors have presented their results for e values up to only 
60 MeV. It would be interesting to see the results of their model also at higher values of e. It 
should, however, be noted that the description of the data for the pn — > put] reaction within 
the three-body model model is less satisfactory in comparison to that for the pp — > ppn 
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FIG. 6: The total cross section for the p + p — > p+p + rj reaction as a function of the excess energy 
with only the pp FSI effects. The solid and dashed curves represent cross section obtained with pp 



FSI calculated with Jost function method and that with the method described in Ref. 



6l| 



reaction. 

We see that the ELM is able to describe both the energy dependence and the absolute 
magnitudes of the experimental cross sections for the r\ meson production in both pp and 
pn channels for excess energies > 15 MeV. However, it underpredicts the pp channel data 
for e values below 15 MeV. Such a trend has also been seen in calculations presented in 
Refs. (29), 0] where the rj meson production in NN collisions has been investigated within 
a relativistic meson exchange model including the initial state interactions and FSI only in 
the NN sub-system. These authors have attributed the near threshold underestimation of 
the experimental total cross section to the non-inclusion of the rjp FSI in their model. 

Since the ingredients of the primary production amplitude of our model have already 
been checked and fixed by calculations done at higher beam energies where FSI effects are 
absent, the underestimation of the pp — > ppq cross section for very low values of e indicates 
that we need to improve the treatment of the FSI effects. Inclusion of the off shell effects in 
the calculations of FSI is one of the likely improvements. The knowledge about the off shell 
nature of the rjN interaction is still very sparse. However, we can use the results presented in 
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FIG. 7: Contributions of iV*(1535) (dashed line), 7V*(1650) (dotted) and JV*(1710) (dashed-dotted 
line) baryonic resonances to the total cross section for pp — ► pprj reaction. Their coherent sum is 
shown by the solid line. 



Ref. [6l| to investigate the effects of using off shell pp FSI on the near threshold rj production 
cross sections. In Fig. 6, we have compared the results for the total cross sections for the 
pp — > pprj reaction obtained by including FSI in pp substate calculated within Jost function 



technique and that obtained with the method described in Ref. |6JJ which includes off-shell 
effects. In the results presented in this figure, no r\p FSI has been considered. We note that 
the off shell effects in the pp FSI do increase the cross section for e < 60 MeV. However, 
it is not enough to explain the underprediction of the experimental data by our theory at 
smaller energies. One needs to have a better understanding of the rjN scattering amplitude. 
Further improvement may come by including the three-body terms in the expansion of the 



scattering amplitude given by Eq. (23). It has been shown in Ref. [32j that NNr/ FSI effects 
calculated within a three-body scattering theory leads to enhanced cross sections at very 
low values of e. 

It should be remarked here that the differences in the cross sections of pn — > pnrj and 
pp — > pprj reactions are not only due to different isospin factors, but also due to differences in 
the FSI effects. The low energy scattering parameters between pp and pn cases are different; 
the latter involves also a triplet spin state together with the singlet one. A crucial difference 
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FIG. 8: Contributions of iV*(1535) (dashed line), iV*(1650) (dotted) and JV*(1710) (dashed/dotted 
line) baryonic resonances to the total cross section for pn — > pnn reaction. Their coherent sum is 
shown by the solid line. 



between them is the Coulomb interaction. This is not included in the three-body model 
calculations of the pp — > ppn reaction reported in Ref . 32| • Inclusion of this term is likely 
to reduce the cross section for beam energies very close to the threshold. 

In Figs 7 and 8 we show the individual contributions of various nucleon resonances to 
the total cross sections of the pp — ► ppr\ and pn — > pnn reactions, respectively, at the near 
threshold beam energies. Similar to the situation at higher beam energies, the cross sections 
are dominated by the iV*(1535) resonance excitation. Since iV*(1535) is the lowest energy 
baryonic resonance having an appreciable branching ratio for the decay into the Nn channel, 
its dominance in this reaction even at beam energies near the r\ production threshold is to be 
expected. The contribution of A?"* (1650) resonance state is small and that of the A r *(1710) 
resonance is even smaller at these lower beam energies. It should, however be noted that the 
resonance-resonance interference terms are not negligible. For the pp — ► pprj reaction, the 
total cross sections are smaller than those corresponding to the A r *(1535) resonance alone. 
For the pn — > pnn reactions the difference between the two is not visible in Fig. 8. 

We found that the inclusion of the amplitudes corresponding to the nucleon intermediate 
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states (the nucleon bremsstrahlung) made a negligible difference in the results reported in 
Figs. 1-8 if the value of the coupling constant gNNr) is taken below 3.0. With the largest 
value of QNNri used in the literature (6.14), the results were affected to the extent of only a 
few percent. This result is in agreement with that reported in Ref. 29]. It is obvious that 
due to a considerable amount of uncertainty in the value of Qnn v (see, e.g. 48, Q, Q); the 



48 



75 



761 ]) and their inclusion 



nucleon excitation amplitudes are quite uncertain (see, e.g. 
makes an insignificant difference to the results reported above. 

In Fig. 9, we show contributions of various meson exchanges to the pp — > ppr] reaction at 
near threshold beam energies. The dashed, dotted, dashed-dotted and dashed-double dotted 
curves represent the contributions of ir, p, a and uo meson exchanges, respectively. Their 
coherent sum is shown by the solid line. We see that the one pion exchange graphs make 
the largest contribution to the reaction in this energy regime. However, a striking feature 
of this figure is that despite a larger value for the Qn*Np coupling used in our calculations, 
the contributions of p meson exchange is still much smaller than that of the pion exchange 



graphs. Hence in contrast to the results reported in Refs. [67J, [73] the p meson exchange 
terms do not dominate the total NNjn production cross sections. To understand this differ- 



ence we note that while in Refs. 



67j, 



771 ] 7 5 7 M couplings have been used for the pNN* vertex, 



we have taken a 750^ coupling for the same which is an extension of the ^NN* couplings 
(due to vector meson dominance reasons 48 Q. This is also compatible with the forms of the 



pNN* couplings used in the literature [46 



78 



791 ] . Since the p meson exchange amplitudes 



calculated with the 750"^ couplings involve delicate cancellations among leading terms, con- 
tributions of this exchange diagrams to the 7] production cross sections are weakened. This 



67 



72]. 



is the main reason for differences between our results and those of Refs. 

In Ref 29J , although the form of the pNN* coupling is the same as ours, relatively lower 



p meson exchange cross sections result from the use of a very small value for the coupling 
constant gN*Np which is based on the lower limit of the branching ratio for the radiative 
decay of this resonance. Our value for this constant, on the other hand, is calculated from 
the branching ratio of the decay of this resonance to Np channel as quoted in Ref. 52 1. 

We see that the u meson exchange process contributes insignificantly to the NNr] produc- 
tion but the u meson exchange terms are relatively more important. Larger contributions 
from the latter has also been seen in other subthreshold reactions analyzed within our model. 
It indicates that a meson exchange may be an efficient means of mediating the large mo- 
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FIG. 9: Contributions of n (dashed line), p (dotted line), u (dashed-double dotted line) and a 
(dashed-dotted line) meson exchange processes to the total cross section for the pp — > ppij reaction 
as a function of the excess energy. Their coherent sum is shown by the solid line. 



80|, 



8l| 



mentum mismatch involve in the meson production in NN collisions 

In fig. 10, we investigate the effects of using PS or PV couplings for the N*Nn(t]) vertices. 
We notice hardly any difference in the cross sections calculated by the two types of couplings. 



Similar results were also observed in Ref 46]. This result is not surprising since the two 



couplings are constructed in such a way that both are equivalent on the mass shell. Of 
course, they start having different energy behavior in the far off shell region where resonance 
contribution is anyway suppressed due to dominance of the corresponding propagator. It is 
only in the NNn case that difference in the PS and PV couplings are obvious with a clear 
preference for the PV coupling in line with the chiral symmetry as discussed earlier. 

After establishing the dynamical content of our model vis a vis the description of the total 
production cross sections, we now turn our attention to more exclusive data. In Fig. 11, 
we show a comparison of our calculations with the data for the angular distribution of rj 
meson in the pp — > ppr] reaction for e values of 15 MeV (upper panel) and 41 MeV (lower 
panel). Since the angular distribution data are normalized to the total cross sections for both 
values of e, we have done the same in our calculations shown in this figure. We note that 
shapes of the angular distributions are described well by our model at both the energies. 
At the lower beam energy, the data as well as our calculations have essentially isotropic 
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FIG. 10: The total cross section for the pp — ► ppr/ reaction calculated with pseudovector (solid line) 
and psuedoscalar (dashed line) couplings for the N*NTr(rj) vertex for the resonances considered in 
this paper, as a function of the excess energy. 



distributions. However, for e = 41 MeV, there is a tendency in our calculations to show 
slight enhancements at forward and backward angles which is typical of the ir exchange 
dominance process in the iV*(1535) excitation. Due to large statistical errors in the data, 
it is difficult to conclude if they show a trend different from our calculations. It will be 
useful to have better quality data with less statistical errors in order to determine if other 
mechanisms which may show a trend different from ours, are also important. In any case, it 
is quite unlikely that the p meson exchange mechanism which might lead to a distribution 



different from that seen in our calculations 
been discussed. 



28j is a dominant mechanism as has already 



IV. SUMMARY AND CONCLUSIONS 



In this paper we investigated the r\ meson production in proton-proton and proton- 
neutron collisions for beam energies ranging from near threshold to about 10 GeV within an 
effective Lagrangian model which has been used previously to describe successfully the pion, 
associated kaon and dilepton production in NN collisions. The interaction between two 
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FIG. 11: Differential cross sections of the pp — » ppr/ reaction as a function of rj meson angle in the 
cm. frame of the total system at the excess energies of 15 MeV (upper panel) and 41 MeV (lower 
panel). The experimental data have been taken from the Ref. [8j. 



nucleons in the initial state is modeled by the effective Lagrangians based on the exchange 
of 7r, p, u and a mesons. The parameters of the corresponding vertices were taken to be 
the same as those used in the previous applications of this model which restricts the free- 
dom of varying the parameters to get fits to the data. The eta meson production proceeds 
via excitation, propagation and decay of iV*(1535), iV*(1650) and iV*(1710) intermediate 
nucleon resonance states. The coupling constants at the resonances-nucleon-meson vertices 
have been determined from the experimental branching ratios of the decays of the reso- 
nances into the relevant channels. Here again we have used the same coupling constants as 
those used in the previous applications of the model at vertices that appeared also in those 
calculations. The interference terms among various amplitudes are included in the total 
transition matrix. 

To describe the data at the near threshold beam energies, the FSI effects among the 
outgoing particles are included by a generalized Watson-Migdal method which allows to 
have these effects in all the three two body sub-systems of the out-going channel. This 
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method involves a parameter which has been determined from the constrained that in the 
limit of FSI in only iViV sub-system the result of the usual Watson method is reproduced. 

In this paper we presented the analysis of the data for total cross sections for the pp — > ppn 
and pn — > pnr] reactions and for the rj angular distributions in the former reaction. With the 
same set of vertex parameters, the model is able to provide a good description of the data 
for the pp — > ppr) reaction at higher as well as near threshold beam energies except for the 
excess energies below 15 MeV where our calculations underpredict the experimental data. 
The experimental total cross sections of the pn — > pnr) reaction are also well described by 
our model. The data for the n angular distributions in the case of pp — > ppn reaction are also 
well reproduced at two beam energies. Imprecise knowledge of the nN scattering amplitudes 
and non-inclusion of the three-body effects are the most likely reasons for underestimation 
of the pp channel data by our model at very low beam energies. 

In this study we have not studied the observables related to meson energy dependence 
e.g., final pp and np invariant mass distributions. There are some open theoretical issues 
concerning the explanation of the corresponding experimental data. While in Ref. 30], 
the inclusion of contributions of the non-s-wave states in the pp subsystem were found to 
be essential to explain these data, the three-body effects in the ppn system and not the 



contribution of 
authors of Ref. 



;he higher partial waves were shown to be crucial for this purpose by the 



32j. At this stage, our theory excludes both these effects. Extension of our 



model to include these mechanisms is vital before we can make some meaningful contribution 
towards settlement of this issue. 

Within our model, one pion exchange processes make the largest contributions to cross 
sections in the entire energy regime. Despite our using a large coupling constant for the 
N*(1535)Np vertex, the cross sections of the p meson exchange process are still lower than 
those of the pion exchange mechanism. Therefore, p meson exchange being the dominant 
mechanism of iV*(1535) resonance excitation 28[ is not supported by our calculations. The 
individual contributions of the u meson exchange diagrams are very small every where. On 
the other hand, the a exchange terms make relatively larger contributions. 

The excitation of the iV*(1535) resonance dominates the NNn production at both higher 
as well as near threshold beam energies. The contributions of iV*(1650) and iV*(1710) 
are small in comparison. However, the interference among various resonance contributions 
is not negligible. Unlike the NNn vertex where there is a clear preference for the PV 
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coupling, the present reaction does not distinguish between PS and PV couplings at the 
N*Nt] vertex involving spin-1/2 even or odd parity resonances. We point out that the 
mechanism of the ppr] production via preferential excitation of the iV*(1535) intermediate 
baryonic resonance state in one-pion-exchange process has received support recently from 
an experimental study [82J of the analyzing powers of the p + p^p + p + rj reaction. 

This work fixes the parameters of the effective Lagrangian model for most of the vertices 
involve in the eta meson production processes. An interesting further check of this model will 
be provided by the analysis of the eta photoproduction data on nucleons (see, e.g. (ssl, 84 ] ) . 
An exciting recent result is that the integrated cross section of the photoproduction of 77 
meson on neutrons shows an additional maximum at center of mass energies around 1.66 
GeV. This has recently been explained in terms of the excitation of the iV*(1650) and 



iV*(1650) resonance states 85(. Furthermore, the vertex parameters derived by us will also 



be useful in applications of effective Lagrangian method in describing the production of 



eta-mesic nuclei in proton and photon induced reactions (see, e.g. [86l. 1871] ). 
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APPENDIX A: FINAL STATE INTERACTION AMPLITUDES FOR THREE- 
PARTICLE STATES 

We give here some clarifications and steps leading to the derivation Eq. (23). 

The total Hamiltonian of the three-particle system is written as H = Hq + U where 
H Q is the kinetic energy operator of the system and the interaction U is taken as U = 
U23 + U31 + U12 = U\ + U 2 + U 3 assuming that the three-particle states interact by means 
of the additive pair interactions represented by Uk- The Green's functions corresponding to 
H and H are, respectively 

G W ( 3 =H ^ __L_ , (Al) 
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G^(E) = lim^o — I (A2) 
hi — Hq ± xe 

We shall also need the Hamiltonian describing the two particles interacting while the third 
one is free, namely, H k = H + U k and the corresponding Green's functions 

The full three-body transition operator T satisfies the Lipmann-Schwinger equations 

T(E) = U + UG T(E), (A4) 

which can also be written as 

T(E) =J2U k + J2 U k G T(E), (A5) 
k k 

Eq.(A5) leads to the iteration 

T(E) = Y,[U k + U k G U k + U k G U k G U k + ...] 

k 

+ £[£4 + U k G U k + ...] Go [Uj + UjGoUj + ...] + ... (A6) 

In the study of the final state interaction problems, one usually has in addition to poten- 
tial U [which is responsible for the transition from the initial state to the free final states 
characterized by the transition matrix To in Eq. (22)], an additional interaction V that 
describes a type of internal interaction among the constituents of the final state. The total 
Hamiltonian is then written as Hq + U + V and one can make use of the standard "two- 



potential formalism" of Goldberger and Watson (see, e.g., |42j) to write the total scattering 
matrix element (Tfi) as as a sum of two terms - one of them involve the matrix elements of 
the interaction U between the exact initial state wave function and the final scattering state 
wave function corresponding to interaction V. From the general theory, we can write 

T fi = < Xf \U + V\4 +) >, (A7) 

where ip^ completely describes the initial state with the outgoing wave boundary condition, 
and Xf is t ne final plane wave state. One can eliminate Xf by introducing wave functions 
0^ ■* which are the eigenfunctions of the Hamiltonian Hq + V as 

= Xf + G^V^ (A8) 
28 



Substituting (A8) into (A7) one gets after some manipulation, 

T fi =< <P ( f ] \U\4 +) > + < ^Mxi >, (A9) 

where Xi is the initial plane wave state. The final state interactions of interest are contained 
in 0j . 

In applications of relevance to us, the second term of Eq. (A9) would vanish because we 
assume V can not create real mesons. Eq. (23) can be obtained by using an iteration similar 
to that given by Eq. (A6) in the remaining (first) term of Eq. (A9). The amplitude To is 
defined in the same way as the first term of Eq. (A9) with a plane wave final state. The 
amplitude Tk in the second term of Eq. (23) represents the matrix elements of the interaction 
Vk between the plane wave and the scattering states of particle ij. We define V& for the 
partition k in the same way as 
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